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Abstract. A 4Mj planet with a 15.8 day orbital period 
has been detected from very precise radial velocity mea- 
surements with the CORALIE echelle spectrograph. A 
second remote and more massive companion has also been 
detected. 

All the planetary companions so far detected in orbit 
closer than 0.08 AU have a parent star with a statisti- 
cally higher metal content compared to the metallicity dis- 
tribution of other stars with planets. Different processes 
occuring during their formation may provide a possible 
explanation for this observation. 
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1. The search for extra-solar planets with 
CORALIE 

Almost 20 planetary companions with minimum 
masses (m2 sin i) less than 10 Mj have so far 
been detected by very high-precision radial veloc- 
ity surveys (Butler et al. 1999 



Marcy et al. 1999a 



Fischer et al. 1999| and Marcy et al. 1999b| for a review 
of older detections). The semi-major axes of their orbits 
range from very small (0.05 AU) to 3AU. Some have 
eccentric orbits, others have secondary more massive com- 
panions, some have both. The large observed spread in 
orbital characteristics of all known planetary candidates 
causes some difficulties in understanding their formation 
process in comparison with our own solar system. It 
also raises the issue of the real nature of these objects, 
particularly the more massive ones. 

In June 1998 we initiated a systematic and large scale 
exoplanet search survey (1600 nearby G and K stars) in 
the southern hemisphere with the new 1.2 m alt-azimuth 
Euler Swiss telescope at La Silla, ESO Chile. The tech- 
nique we are using to detect planets is to look for a stellar 
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reflex motion due to an orbiting planet by very precise ra- 
dial velocity measurements. The CORALIE echelle spec- 
trograph is used to measure star spectra from which the 
Doppler effect is then computed. 

CORALIE is an improved version of the ELODIE 
spectrograph (Barannc et al. 1996) with which, 4 years 



ago, the first extra-solar planet orbiting a star (51 Peg) 
was discovered ( [Mayor fc Queloz 1995D . The CORALIE 
front-end adaptator is located at the Nasmyth focus of 
the Euler telescope. Two sets of two fibers can alterna- 
tively feed the spectrograph which is located in an iso- 
lated and temperature controled room. The set of fibers 
used for high precision radial velocity measurements in- 
cludes a do uble scrambler devi ce designed by Dominique 
Kohler (see [Queloz et al. 19991 for references) to improve 
the stability of the input illumination of the spectrograph. 
Thanks to a slightly different optical combination at the 
entrance of the spectrograph and the use of a 2k by 2k 
CCD camera with smaller pixels (15/im), CORALIE has 
a larger resolution than ELODIE. A resolving power of 
50,000 (A/AA) is observed with a 3 pixel sampling. As 
with the ELODIE spectrograph, CORALIE makes use of 
on-line reduction software that computes the radial veloc- 
ity of stars several minutes after their observation. (See 
Baranne et al. 1996 for details about the reduction pro- 



cess). The simultaneous thorium technique is used to cor- 
rect any instrumental drifts occuring during the star ex- 



posure (see Queloz et al. 1999 for details). The many im- 
provements carried out in the thermal control and the 
resolution of the instrument, as well as in the reduction 
software, yield a factor two improvement in the instrument 
precision compared with ELODIE. 

2. A planet orbiting Gliese86 

Gliese86 (HD13445, HIC 10138) is a bright {my = 
6.12) early K dwarf (B ~ V = 81, T^q = 5350 K, 
\og{g) = 4.6, [Flynn fc Morell 1997 ) from the southern 
hemisphere, in the Eridanus (River) constellation. It is 
a close star, 10.9 pc away from our Sun (tt — 91.6 mas, 
measured by the Hipparcos satellite). Its absolute magni- 
tude is 6.257, yielding (with BC = —0.2) a luminosity 
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Table 1. Orbital elements of Gliese 86 after correction of 
the 0.36 ms~^ d~^linear drift of the 7-point. 
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L — 0.27 Lq. It is a high proper motion star, slightly 
metal poor ([Fe/H]= -0.24, |Flynn fc MorcU 1997[ ). It 



has low chromospheric activity (logi?'jjj^ ~ —4.74, 
Osten 1997[). No rota tional broadening has been 



Saar 



detected ( ^aar fc Osten 1997 ) and there is only an upper 
limit on the Li content in its atmosphere (N(Li)< —0.24, 
Favata et al. 1997). From Hipparcos photometry, the star 
is stable (cr(IIp) — 0.008). In summary, GHese86 bears all 
the characteristics of a few billion year old K dwarf from 
the old disk population. In the H-R diagram, Gliese 86 lies 
slightly below the ZAMS. However we believe that there 
are enough uncertainties in the temperature and bolomet- 
ric correction estimates of Gliese 86 - stemming from its 
low metal content - to believe that its location in the HR 
diagram below the ZAMS is not significant. 

A 15.8 day period radial velocity variation has been 
detected from CORALIE measurements (Fig. 1). In Ta- 
ble 1 are listed the orbital elements of the best fit so- 
lution (least square) for an orbital motion after correc- 
tion of a 0.36 ms~^ d^^linear drift (see below). Assuming 
a 0.8Mq for the primary and that the radial velocity 
effect is caused by the orbital motion of the star, we con- 
clude that a 4 M j companion (minimum mass) is orbiting 
Gliese 86. 

The planetary companion to Gliese 86 is close to its 
host star with a 0.11 AU semi-major orbital axis. It has 
a low, although 99% significant non-zero, eccentricity 
(Lucy & Sweeney 1971). The 7ms~^residual from the fit 



indicates very low intrinsic intrumental errors from night 
to night, taking into account that each measurement has 
approximately 5 ms~^ photon noise error and could be af- 
fected as well by some low level radial velocity variations 
intrinsic to the stellar atmosphere. Such low instrumen- 
tal error agrees with the instrumental error measured by 
P(x^) analysis of all the stars of our samp le so far observed 



(about 300). See Duquennoy et al. 1991 



for a detailed de- 
scription of the instrument error estimate by the P{x^) 
statistic. 

A long term drift of the radial velocity (0.5 
ms~^d~^) is observed from 20 years of CORAVEL 
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Fig. 1. Phased orbital motion of Gliese 86 corrected from 
the long term drift. The solid line is the best fit orbit. See 
orbital elements in Table 1 



measurements (Fig. 2). With the 300ms^^typical pre- 
cision of CORAVEL radial velocities, the short orbit is 
marginally detected in the last measurements. Interest- 
ingly, with the recent CORALIE measurements a smaller 
0.36 ms~^ d~^drift is observed (Fig. 3). A statistical anal- 
ysis of the reliability of the drift correction shows that 
an orbital solution without drift correction has 0.0001% 
chance to occur (x^ k, 210). The probability jumps to 40% 
when the linear drift correction is taken into account. A 
conservative 7 ms^^ instrumental error is assumed for this 
calculation. 

The period measurement of the short period planetary 
companion is still not accurate enough to correct the old 
CORAVEL data from their extra scattering and obtain 
a precise drift estimate from these measurements. Thus, 
the difference in drift slope between the old CORAVEL 
measurements and the recent CORALIE measurements is 
perhaps significant, but remains to be confirmed by fur- 
ther measurements during the course of the next season. 

The long term radial velocity variation is the signature 
of a remote and more massive companion. The use of his- 
torical radial velocity data together with the CORAVEL 
and the CORALIE observed drifts suggest a stellar com- 
panion with a period longer than 100 yr (semi-major axis 
larger than 20 AU). A direct detection would be worth 
attempting since the star is close to us. 

Alternative explanations to a low mass companion 
to explain the observed 15.8 day period radial velocity 
chan ge of Gliese 86 would be activity related phenom- 
ena ( Saar fc Donahue 1997 ). However Gliese 86 doesn't 
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Fig. 2. Filled dots: radial velocity drift observed with 
CORAVEL. A mean 0.5 m s^^ d^^ variation of the radial 
velocity of Gliese86is measured {dotted line). The 15.8 day 
reflex motion from the planet is marginaly seen as an ex- 
tra scattering in the last CORAVEL measurements. In the 
right lower box is displayed for comparison the 7-point 
drift measured with CORALIE. (Note that the time scale 
is artificially extended for the sake of a better display). 
Open dots: previous measurements found in the literature. 
No error bars are displayed but a typical 2-5kms~^ error 
may be assumed for these measurements 
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exhibit any of the classical activity signatures seen on 
young stars, as for exemple HD166435 (Queloz et al. in 
prep). Gliese86 has no chromospheric activity. No rota- 
tional broadening is detected either, and its photometry 
is very stable. Therefore, the planetary hypothesis is most 
likely the correct interpretation for the observed periodic 
radial velocity changes. 



Fig. 3. Top: CORALIE observed 7-point drift (residu- 
als from the short orbit fit). The dashed line is the best 
linear fit after correction of the short period from the 
planetary companion. It corresponds to a linear drift of 
0.36ms~^d~^. The dotted line is the 0.5ms^^ d~^drift 
measured over 20 yr CORAVEL data. Bottom: Set of 
3 non-variable stars measured at the same period of 
time (open dots: HD10700, stars: HD39091, filled dots: 
HD67199. No instrument zero point drift is observed 



3. Discussion 

The observed orbital charateristics of planets are the di- 
rect outcome of their formation processes and of their 
evolutions. Therefore, these characteristics may be used 
to retrace their formation mechanisms and to constrain 
theories of planetary formation. The recent spectroscopic 
studies of stars where planets have been detected have 
shown that the host star itself may also bear marks 
from some processes occuring during planetary formation 
( ponzalcz 1997 ). More specifically a large number of plan- 
ets with short-period orbits have surprinsingly metal rich 
host stars. These planets, very close to their stars, are usu- 
ally referred as 51 Peg like, or "hot Jupiters". The metal- 
licity of their host star is much higher than the "average 
field star" and is no t the result of a sel ection process in 
the survey samples ( Marcy et al. 19991 ). Typical metal- 
licities similar to field stars may be assumed for the stars 



from various surveys, since these star samples have not 
been selected from any metallicity criteria. 

If we look in more detail at all the planets with semi- 
major axis less than lAU, where the number of detec- 
tions is significant and not strongly high-mass biased, we 
observe a relation between the semi-major axes of the 
planetary orbits and the metal content of their host stars. 
All planets with semi-major axes less than 0.08 AU seem 
to have a star with an unusually very high metal con- 
tent compared to other stars with planets (sec Fig. 4) . 
Actually, a comparison of the two distribution using the 
Kolmogorov-Smirnov test indicates a 99% probability that 
the two distributions are indeed different. 

The unusual metal content of the short orbit planets 
had been pointed out sh ortly after the detection of 51 Peg 
( Mayor fc Queloz 1995 ). But now, with the large number 
of detections of similar systems and others with slightly 
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Fig. 4. Top [Fe/H] content of stars versus semi-major axis 
of the planetary orbit for all known planetary candidates 
with a 7712 sin i < 12 Mj and a < lAU. [Fe/H] measure- 
ments are from Gonzalez 1998 and [Gonzalez et al. 1999 . 
A typical 0.06 error are given by the authors. The dot- 
ted line connects the dots representatives of the two inner 
planets orbiting the star v Andromedae. 
Bottom Distribution of the metallicity of stars with a 
planet (solid line) . The shaded area indicates stars with a 
planet closer than 0.08 AU. Note the wAnd multiple sys- 
tem and the planet orbiting the late M star GJ876 are not 
included in the histogram (open dots on Top diagram) 



larger semi-major axes, we observe a typical distance (or 
period) for which this unusual high metal content is sys- 
tematically observed. A possible explanation may be re- 
lated to some very specific processes occuring during the 
formation of these very close systems. However the large 
incertainties on the estimation of the age of these systems 
and the small mass range of primaries are noteworthy. 
Therefore it is difficult to completely rule out a stellar 
population effect. 
The 



migration theory (Lin & Papaloizou 1986, 
Lin et al. 1996| [Ward 1997] , [Trilling et al. 1998D is one of 
theories that has been called for to explain the existence 
of very close planets that were not described b y the "clas- 
sic" solar-system planetary formation model (Boss 1995, 



Lissauer 1995). But so far, we have a poor understanding 
of the way the planet stops its migration. The two 
different metallicity distributions pointed out in this 
article are perhaps a new clue to a better understanding 
of the migration process or the likelyhood of an in-situ 
formation ( Bodenhcimer et al. 99| ). 

Others scenarios involving strong gravity inter- 
actions with other planets have been proposed 
as a possible origin for small planetary orbits 



(WeidenschiUing fc Marzari 1996, Rasio et al. 1996) 



Since these models are purely driven by dynamical inter- 
actions it seems a priori difficult to expect any metallicity 
enhancement effect. Moreover such scenarios do not 
really explain the very small orbit planets like 51 Peg. 
However, if one believes that the high metal content of 
the s tar is the end-result of a planet swallowed by the 
star ( Sandquist et al. 1998 ), the gravity interaction is a 
possible means to send planets into their stars. 

The precision of surveys from which the planets have 
been found so far has been limited to the detection of 
systems with a yr-amplitude (K) larger than 25ms~^. 
Therefore it is still premature to compare the mass of the 
two sets of planets because there is a direct relationship 
between the amplitude of the radial velocity curve, the 
semi- major axis, and the minimum mass of the planet that 
can be detected. However if we limit our comparison to a 
sample free of such bias, including planets only having 
semi-major axes smaller than 0.3 AU, we may be inclined 
to believe that metal rich stars tend to host on average 
less massive planets than solar stars or metal poor stars. 

The distribution of mass of all planets that have been 
detected can also be studied with a restricted sample of 
planetary systems in order to avoid a biased selection to- 
wards small orbits and massive systems. We see that with 
a sample restricted to minimum masses greater than 1 Mj 
and semi-major axes smaller than 1.3 AU, the number of 
planets per mass bin is almost constant from 1 to 5 Mj 
and then drops suddenly for more massive companions. 
This reinforces the idea that a maximum planet m.ass may 
lie somewhere close to 5-7 Mj as pointed out earlier by 
Mayor et al. 9^ 

New discoveries and improved detection precision will 
allow us to get a better picture of the relation between 
the mass and certain orbital characteritics of planets and 
some peculiarities seen in the atmosphere of their host 
stars. This will perhaps enhance our understanding of the 
mechanisms of planetary formation. 
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